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Abstract 

The effective action on an orbifolded sphere is computed for minimally coupled 
scalar fields. The results are presented in terms of derivatives of Barnes ^-functions 
and it is shown how these may be evaluated. Numerical values are shown. An 
analytical, heat-kernel derivation of the Cesaro-Fedorov formula for the number of 
symmetry planes of a regular solid is also presented. 
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1. Introduction. 

In earlier work [1] we have shown that the ^-function, Cr{s), on orbifold- factored 
spheres, S'^/F, for a conformally coupled scalar field, is given by a Barnes C^function, 
[2], Cd{s, a I d), where the di are the degrees associated with the tiling group F. The 
free-field Casimir energy on the space-time IRxS'^/F was given as the value of the C- 
function at a negativ e integer which evaluated to a generalised Bernoulli function. 
In the present work we wish to consider the effective action on orbifolds S^/F which 
this time are to be looked upon as Euclidean space-times. In particular we will 
discuss d — 2 and d — 3, concentrating on the former. 

The simplifying assumption in our previous work was that of conformal cou- 
pling on IRxS'^/F. This made the relevant eigenvalues perfect squares and allowed 
us to use known generating functions to incorporate the degeneracies. From the 
point of view of field theories on the space-times S'^/F, retaining this assumption 
would be rather artificial. A more appropriate choice would be minimal coupling, 
or possibly conformal coupling, on S^/T. (These coincide for d = 2.) 

The quantities in which we are interested are Cf(0) Cr(0). The latter 
determines the divergence in the effective action and the former is, up to a factor 
and a finite addition, the renormalised effective action [i.e. half the logarithm of 
the functional determinant). 

2. Eigenvalues, degeneracies and zeta functions. 

For the aforementioned conformal coupling, the eigenvalues of the second order 
operator -A2 + = {d - l)/4d) are 



with degeneracies that we shall leave unspecified here. 

In our previous work [1] we showed that the corresponding Neumann and 
Dirichlet ^-functions on S^/T were. 




(1) 



Cr(^) = Cd(2.,(rf-l)/2|d), 
Cf^(«) = Cd(2s,Ec^i-(c^-l)/2|d), 



(2) 



(3) 



where the general definition of the Barnes ^-function is 




(4) 



This shows that the eigenvalues are given specifically by 



An = (o + Hl.d)^ 



(5) 
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the degeneracies coming from coincidences. The parameter a is (cZ — l)/2 in the 
Neumann case and comparison with the previous form shows that the integer n — 
2m.d + 1, m = upwards. For Dirichlet conditions, a = ^ c/^ — (d — l)/2 and then 
n = 2m.d — 1 with m = (1, 1) upwards. The interpretation in two dimensions is 
that the angular momentum is L = m.d for Neumann and m.d — 1 for Dirichlet 
conditions. 

Turning to minimal coupling, = 0), the eigenvalues of the Laplacian are 

A, = (a + m.d)2-i^^. (6) 
and the corresponding (^-function is 

C(s) = y"7 ra- (7) 

The origin m = is to be omitted for Neumann conditions, when the ^-function is 
denoted by C{s). 

Consider a sum of the form 

^ [{a + m.d)^ - a^) 

so that 

C{s) = as)-{a^-a^)-'. (9) 

For minimal coupling, a = (d — l)/2, while for conformal coupling in d-dimensions, 
a = 1/2. We concentrate on minimal coupling. 

A standard way of obtaining information about an expression such as (8) is 
to perform a binomial expansion to produce a sum of known ^-functions, in the 
present case a sum of Barnes C-functions, 

= g a^r s{s + l)...{s + r-l) ^^^^^ ^ ^ I ^^^^ 

r=0 ^' 

From this, the value of ({s) at a nonpositive integer is easily found. For example 
the important value ((0) is given by 

C(0) = Cd(0, a I d) + - ^ —N,r{d) (11) 

2 r 

r=l 

and, more generally, we have 
a-n) = ±{-ay (!) U2r - 2n, a\d)+^-=^± !i^!_zil!„^n+-jV.., (12) 



where u = d/2 if d is even and u= (d — l)/2 if d is odd. 
Nr{d) is the residue defined by 

Cd(s + r,a\d)^^^^ + Rr(d) as s ^ 0, (13) 
s 

where 1 < r < d. Expressions for the residue and remainder involve generahsed 
Bernoulli functions and can be found in Barnes [2] . For shortness, their dependence 
on the parameter a is not indicated. 

The form of the residues given by Barnes [2] is 

where dS^^^^ (o) is the (r + l)-th derivative of Barnes' generalised Bernoulli poly- 
nomial d'S'i(a). The general relation with the more usual polynomials, [5], will not 
be given here. Specific forms are 

Barnes also gives the values 

U-n, a I d) = .S^Ua) = (jg^B«t(a | d). (15) 

From (11), (14) and (15) we find, for two dimensions, 

C (0) = (0) = (3 - 3((ii + d2) + (dl + d2f + did2) . (16) 

This corrects our previous expression [1] obtained by an incorrect manipulation 
of the heat-kernel. 

3. The derivative of the zeta function. 

The derivative at s = is a little more difficult to find. From (10) a first step is 
C{0) = 2CA0,a\d) + J2^[R2r + ^N2rJ2l]+ E ^Cd(2r,a|d). (17) 

r=l \ 1 / r=u+l 
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The integral representation of the Barnes ^-function allows the final sum in 
(17) to be written as 



V r dr ^ ^^^p(-^^) 

2 / exp(-aT) coshar- > \ ^ , ^ -v- (18) 

k \ ^0 ^2r)! ^ rn.(l-cxp(-d,r)) ^ ^ 

In the Neumann case a = a and there is an infra-red, logarithmic divergence at 
infinity caused by the zero mode which will be taken care of by the transition to 
(9). 

Although the integral converges nicely at r = 0, the individual terms of the in- 
tegrand do not. It is enough to introduce another ultra-violet analytic regularisation 
and define the intermediate quantity, 



2 / exp(-aT) coshar- > \ ——, (19) 

^ ^ (2r)! ; n. (1 - exp(-rf^r)) ^ ^ 



whose s = limit gives (18). 

After continuation, (19) integrates to 



Vis) (Cd(s, a-a\d) + Cd(s, a + a | d)) - 2 J] ^r(s + 2r)Cd(s + 2r, a | d). (20) 



r=0 



As s tends to zero, each term in (20) yields a pole and a finite remainder. The 
poles must cancel and so 



Cd(0, a-a\d) + Cd(0, a + a I d) - 2Cd(0, a | d) = ^ —N2rid). (21) 

This condition is an identity between generalised Bernoulli functions. Combinined 
with (11) it produces the symmetrical expression 

C(0) = I (Cd(0, a-a\d) + Cd(0, a + a\d)). (22) 
The finite remainder in (20) is 

CM a-a\d) + CM a + a\d)- 2(^(0, a\d)-J2 —R2M- 



r 
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u 2r 

7(Cd(0, a-a\d) + Cd(0, a + a\d)- 2Cd(0, a | d)) - ^ —i;{2r)N2r{d) (23) 

r=l 

which, in view of (21), can be written 

C^(0, a-a\d) + C^(0, a + a I d) - 2(^(0, a | d)- 

J] ^i?2r(ci) ^ 5^ ^(V(2r) + 7)iV2.(c/). (24) 
Combining this with (17) we have finally 

u 2r 

C'(0) = C^(0,a-a I d)+C^(0,a+a | d)-^ ^ (2V^(2r)-V'(r-l)+7)iV2r(ci). (25) 

The fact that the remainders have cancelled, suggests that there is a more rapid 
route to this result. 

Apart from the final term, (25) is the expression that would have been obtained 
by a naive application of the 'surrogate' ^-function method which is based on the 
product nature of the eigenvalues, {a — a + m.d) {a + a + m.d), in (8) followed by 
an application of the rule Indet (AB) = Indet A + lndet B. This method is suspect, 
as discussed by Allen [3] and by Chodos and Myers [4] . Allen [3] derives a particular 
'correction' term as in (25). He also points out that (22) could be expected on the 
basis of the eigenvalue factorisation, being the average of the regularised dimensions 
of the operator factors. 

The final term in (25) can be rewritten 



E ^ (2^(20 - ^(r - 1) + l)N2r{d) = X: ^N^rid) J2 

r=l r—1 

In order to evaluate the effective action we must substitute the appropriate 
values of a and a for Neumann and Dirichlet conditions into (25). In the former 
case it is also necessary to remove the zero mode i. e. to use i^. The relevant quantity 
then is the F-modular form p, defined by, [2] , 

limC(0,e I d) = -lne-lnp^(d). (26) 

e— i-O 

We find the following basic expressions 

u 2r r — 1 ^ 

C(0) = - Inpd(d) + C^(0, d-l\d)- ln{d - 1) - E ^N2r{d) ^ (27) 
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and 



u 2r 



C;(0) = C^(0, do\d) + CM do + d-l\d)-J2 ^N,r{d) J2 ^j-rj^ (28) 

r=l 

where do = di — d + 1 is the number of reflecting planes in F. We recall that 
a = {d—l)/2 for minimal coupling and that in (27), N is evaluated at a = {d — l)/2 
and in (28) at a = do + {d - l)/2. 

In the case of the two-sphere, (27) and (28) become 

C;(0) = -lnp2(d)+C^(0,l|d)-^ (29) 

and 

C(0) = C2(0, do\d)+ C^(0, do + 1 I d) - ^ (30) 

where we have set g = did2, the order of the rotational part of T. 
For the three-sphere 

C;(0) = - lnp3(d) + C3(0, 2 I d) - In 2 + 1^ (31) 

and 

(0) - C3(0, do I d) + C^(0, do + 2\d)-^ (32) 

where g = did2ds. We note the change of sign in the last term. 

Equations (29) to (32) are the calculational formulae we shall use in the rest of 
this paper. It is also possible to evaluate the derivative of the C"function at negative 
integers, ('{—n). This would be relevant if we were interested in the effective action 
on product spaces like R x R'^ x S'^/F. A few details are presented in the appendix. 

Although our main interest is in minimal coupling, it should be mentioned 
that the result (25) can be used immediately for massive fields, assuming that the 
appropriate value of a is real. This means that the mass k is restricted to the region 
< K < {d ~ l)/2. For larger masses a slightly different continuation is needed. 

4. The derivative of the Barnes zeta function. 

We turn now to the evaluation of the derivatives needed in (29) and (30). A pre- 
liminary step is to remove any common factors of the degrees di and d2 by setting 
di = CBi with ei and 62 coprime so that the denominator function in (4) equals 
c{h + m.e) where h = a/c. 

The summation in (4) is rewitten by introducing the residue classes with respect 
to e. On setting 

mi = 77-162 TO2 = n2ei -l-pi. 
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where < < — 1, the denominator function in (4) equals 0(6 + 6162(^1 + 712) + 
P2ei +P162) and the sum over m becomes 



00 ^ 
X - \ ^ 1 + n 



C2(s, a|d) = c > > -jz ^ 

V> + ei62n + p2ei + P162)* 



p n=0 

^ E (^T^ + (^) Ed - E ^d^. 



(33) 



where 

A' = 6ie2n +^261 +^162, and = 1 1 . 

6162 61 62 

Consider the integer N = eie2n + P2ei +P162. As n ranges over to 00, 
and the pi over their domains, will likewise run over this infinite range with 
the exception of some integers < 6162 at the beginning, specifically those integers 
that equal p2ei + pie2 mod 6162 for ^2^1 +Pie2 > 6162- We denote these missing 
integers by Vi. Apart from these terms, the first sum in the second line of (33) will 
immediately give a single Hurwitz ^-function, 

C2(., a I d) = ^ (^Cr{s -l,b)-Yl (b + liY-^ ) + (^) ~ 

(34) 

(34) is a convenient form for numerical evaluation. It provides an explicit 
analytical continuation of this integral Barnes (^-function. 

For the derivative at s = we find, after inserting the known values of the 
Hurwitz ^-function and its derivative, 

C2 



C^(0, a I d) = - \^C'n{-l, b) + Y,{v, + h) + h) 

- J (^CR{-l,b) - + Inc (35) 

+ 5^(1 - Wb) (in (r(«;b)/^/(27r)) - (1/2 - w^) Hd/c)) • 
p 

Letting a tend to zero in (35) and comparing with the definition of Inp, (26), 
one finds that 

lnp2(d) = -y (^C'ji{-1) + J2^i'^^^i^ -7 + 

- 5]'(1 - wo) [in (r(«;o)/ vTStt)) - (^ - wo) ln(^/c)) - ^ ln(^/27rc) 
p ^ ^ 

(36) 



where wq = Pi/ei ■A-'p^l^i and the dash means that the term pi = p2 = is to be 
omitted from the sum. 

Barnes gives a formula in terms of the multiple F-function, 



(38) 



Formal expressions for the functional determinants are thus 

r.(i) 

and 

p-Cl, (0) _ ^2(d) .ggX 

F2(do)F2(do + l)' ^ ^ 

Our results, (34) and (35), (36), can be thought of as computational 
formulae for these functions in terms of simpler ones. 

It is not necessary to rearrange the summation as in (33). We have done so in 
order to extract the term C/j(s — 1, 6). If the summation is left as in the first line of 
(33), it can immediately be turned into a sum of Hurwitz ^-functions, 



C2(s,a I d) = ^ (Cfl(s - IjW^b) + (1 - w^6)Cfl(s>«^6))- 



(40) 



p 

Then we have the alternative form 



(41) 



C^(0, a I d) = ln(c/^) (Ci?(-1, -^h) + (1 - «'6)Ci?(0, w^)) 

p 

p 

= (6a2 - 6a(do + 1) + (c?o + 1)' + 9) M^ld) 

+ E ^6) + (1 - w^) In (F(«;5)/y(27r))) . 

p 

In this way we do not need to find the missing integers (nor even the common factor 
c) but the price is the multiple evaluation of C^( — 1, Wh) by a numerical procedure. 
There is no difficulty in this but (34) is faster and more accurate. Equation 
(41) constitutes a useful check. 

5. The point groups. 

A limited test of our formulae is provided by the dihedral case, F = \q\ in Coxeter's 
notation [7,8]. (Schonflies would write Cqy and it is Cg[Dq in Polya and Meyer 
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[9,10]. Table 2 in [11] has a complete list of equivalents). The degrees are di = q, 
d2 = 1, so c = 1, g = q and do = q. There are no missing integers i^j and, 
furthermore, p2 = 0. The fundamental domain is the lune, or digon, (qql). 

For q — 1 there is a single, equatorial reflection plane, the fundamental domain 
being a hemisphere, (111) (a spherical triangle with every angle equal to tt). An 
alternative notation for this domain is Ai, [7]. In this extreme case, pi is also zero 
and the expressions rapidly collapse to 

lnp2(l, 1) = -Ck(-l) - In ^27r), (2(0, 1 | 1, 1) = Cfl(-l) 

and 

C'(0,2|l,l)=Ci^(-l)+ln^/(27r). 
Thus, on the hemisphere, from (29) and (30), 

^(0) = 2C^{-1) - In ^{2n) - i, C^(0) = 2C^(-1) + In ^{2n) - \, (42) 

which agree with the results exhibited by Weisberger [12]. Our value of CAr(O) = 1/6 
does not agree with [12]. 

The sum of the Neumann and Dirichlet expressions should reduce to the full- 
sphere result derived by e.g. Hortagsu, Rothe and Schroer [13] and later by Weis- 
berger [14] . We find 

Cs^(O) = 4C^(-1) - ^ ~ -1.161684575 

agreeing with these earlier calculations. There are many discussions on spheres 
bounded equatorially by spheres. 

We give the explicit formulae for the next value of q, q = 2, corresponding to 
a quarter sphere, 

Cn{0) = Cr{-1) - In ^/(27^) - i, C^(0) = (^^(-1) + In ^/(27^) - ^. (43) 

Adding these expressions gives half the full-sphere value. 

The results for higher values of q are shown in Fig.l, where we plot the effective 
action W = — C'(0)/2. It is shown in the appendix that 

G(0)-C^(0) = -ln(27r) (44) 

for all [q], as born out by the numbers. 

For completeness wc record the values of C(O) obtained from (16), relevant for 
the conformal anomaly, 

C(0) = C.(0) = ^(l + 9^). (45) 
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We turn now to the extended dihedral group, [q, 2] , of order 4q', obtained from 
[q] by adding a perpendicular reflection. It is the complete symmetry group of the 
dihedron. (In [9,10] this group is D^j {q even) and Dq[D2q {q odd). The Sch5nflies 
equivalent is Dqd-) 

If q is odd, c = 1, di = q, d2 = 2, do = q + 1 and g — 2q, while, if q is even, 
c = 2, ei = q/2 and 62 = 1. For odd q, the missing integers are 1, 3, . . . , — 2. 
There are no missing integers if q is even. 

The fundamental domain is the spherical triangle (22^). When q = 1 this 
domain is the quartersphere lune and the results coincide with (43). The group 
isomorphisms are [1, 2] ^ [2] ^ [1] x [1] (or CafDa = DifDa since C2 = Di). 

Generally one has [q,2] ^ [q] x [1], in particular, [2,2] ^ [1] x [1] x [1] which 
corresponds to three perpendicular reflections with the eighthsphere, (222), as fun- 
damental domain. 

The hemisphere, quartersphere and eighthsphere are the intersections of 
with (R+xM^), (M+xlR+xlR) and (R+xR+x]R+), respectively. The positive real 
axis, R"*", is the positive root space of the SU(2) algebra, Ai (c/ [6]). Fig.2 displays 
values of W for bigger orders. 

The rotation part of [q, 2] is the complete symmetry group of the regular q-gon, 
{q}, and is the dihedral group in its guise as a group of rotations. Coxeter denotes 
it by [g, 2]+ and Polya and Meyer by Dq. As stated, its structure is [g, 2]+ = Dq. 
When q is odd there is the curious isomorphism [2q] = [2,q]. 

It is only a matter of substitution to work out the the values of (29) and (30) 
for the other reflection groups which are the complete symmetry groups of the 
spherical tessellations {3,3}, {3,4} and {3,5}. We flnd — C'(0)/2 for (Dirichlet, 
Neumann)-conditions to be (0.45603,-0.34216) for = [3,3], (0.2508,0.001915) 
for Oh = [3, 4] and (-0.10538, 0.45014) for h = [3, 5]. 

The fundamental domain of [p, q] is the spherical triangle (pqr). The rotational 
part of [p, q], i.e. [p,q]~^, is often denoted by {p, q,r). 

7. The Cesaro-Fedorov formula. 

It is interesting to check the formula (16) by remembering that C(0) is a local object 
related to the constant term in the short-time expansion of the heat-kernel. The 
general formula for a two-dimensional domain, J\4, with boundary dM. = [j dAii 
is 

= i / + / ^(O^^ + ^E^^^ (46) 

247rJ^ 12^Jqm. 247r^ a 

where the a sum runs over all inward facing angles at the corners of dA4. 

In the present case R = 2 and the extrinsic curvature, k, vanishes since the 
boundaries of the fundamental domains are geodesic. Therefore 
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where we have used the standard formula for the area of a spherical triangle. This 
agrees with (16) if the formula 

2dQ{do -l)^g{p + q + r-3) (48) 

is taken into account. In fact our derivation can be thought of as an analytical proof 
of this relation which is a slight generalisation of equation 4-51 in [7]. (Coxeter has 
r = 2 and g = 2Ni.) 

Coxeter indicates a purely geometric proof and points out that (48) is equivalent 
to a formula discovered numerologically by Cesaro [15] and is a special case of an 
earlier result of Fedorov [see 16]. (48) is virtually identical to the equation on pl77 
of [15] with the correspondances X = do, n = r, p — p and q = q- An extension 
to higher dimensions is possible using the generalisation of (46) that includes the 
results of Fedosov on polyhedral domains, [17]. 

7. Scaling and limits. 

The results given so far are for a unit sphere. For radius R, simple scaling gives the 
relation 

C'(0;i?)=C'(0)+21ni?C(0) (49) 

where C'(0) = C'(0; 1) and C(0) = C(0; 1) = C(0;i?). 

The effective action should incorporate an arbitrary scaling length, L, by 

Wl = -^C'(0; R) + InL C(0) = -^C'(O) + HL/R) C(0). 

The figures show just the first term. 

Consider the dihedral case [q] and let q and R tend to infinity in such a way that 
the equatorial width of the fundamental domain (qql) remains fixed at = irR/q. 
From (49) and (45), whence C(O) ?/12, we have 

C'(0;i?)^ limC'(0) + ^ln(^). (50) 

The area of {qql) is Aq — Ifp^qj-K and requiring the density, C'(0;-R)/^q, to 
remain finite as — > oo entails the leading behaviour 

C'(0)^-^lng + O(g). (51) 

Numerically we find 

C'(0) ^ -|lng + 0.497509gf« -|ln(g/19.79) (52) 
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so that the density becomes 



^^->Y^ln(6.299/?). (53) 

Geometrically, it might be imagined that in the limit R — oo, since the sphere 
becomes flat, the rescaled Imie, (ooool), would be an infinite strip of width /3. 
Defining the strip coordinates x — Rcf) and y — R{n/2 — 9), the spherical Laplacian 
does become the usual Cartesian one as — > oo. However, the influence of the 
infinitely sharp corners at the poles persists, even though they are infinitely distant, 
producing an anomaly density of 7r/12/3^. On the rectangular strip, infinite or not, 
the integrated anomaly equals 1/2 and so the density vanishes in the infinite case. 

8. The three-sphere. 

The expressions for the three-sphere are (31) and (32). Then we require, 

CsiSja I d) = 7 ; ; 

(a + midi + 777,202 + mads)* 

We will not attempt to extract a single ^-function as we did previously but will just 
reduce the sum to a finite one over Hurwitz ^-functions in a not very symmetrical 
nor economic fashion. 
The residue classes 

1712 = din2+pi, mi=d2ni+p2 

are introduced so that the denominator function reads (a + did2{ni + 772) +P2<^i + 
Pid2 + m3ds). The sums over ?7i and ?72 can be transformed by defining n = ni +772 
and doing the sum over ni — ?72 to yield the intermediate form 

00 ^ 

C (s a\ d) = 

2^ 2^ (a + did2n+p2di+pid2 + m3d3y' 

Pi,P2 n,m3=0 ^ f f o J/ 

The further residue classes 

ms = dxd2nz + ps, n = ^3774 + p\ 
are introduced and the sum and difference defined by 

77+ = 774 + 7l3, 7l_ = 7l4 — Uz. 

The denominator is independent of 77_ while the numerator equals 1 + d^irij^ + 
n-)/2 +]34. Since the range of 77_ is symmetrical about zero (from — ri_|_ to 77_|_ in 
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steps of 2) the n_ term gives nothing and there is a factor of (1 + n+) multiplying 
the rest. The sum may therefore be written 

,i \ A\ V- (l + ri)(H-c^3n/2+p4) 

C3(^.«|d)-E ^f^. 

where p = (pi,P2,P3,P4), / = o + dxd2V\ + Vi^x + Pid2 + Psds, g = did2dz and we 
have set n = ns^ for notational simplicity. 
The numerator is reorganised to 

(1 + n)(l + ^3^/2 + = A (F + G{f + gn) + (/ + gnf) 

^9 

where 

F = {A-g + did2P4){A - did2P4 - 2did2), G = g + 2did2 - 2A 

with A being the combination A = a + diP2 + d2Pi + (isPs- 

Thus, finally, we arrive at a finite sum of Hurwitz ^-functions, 

C3(., a\d) = ^Y. [4 -) + -^M' - 1, -) + 4=2Cr{s - 2, ^) 1 (54) 

2^ ^ i9 9 g g g g \ 

which constitutes a possible, but inefficient, continuation of the Barnes ^-function. 



9. The honeycomb groups. 

The three-dimensional analogues of the polyhedral tessellations, {p, q], of the two- 
sphere are the spherical honeycombs {p, g, r}, [7,8,18]. The reflection groups [p, g, r] 
are their complete symmetry groups, the fundamental domains being subspaces of 
the honeycomb cells. A numerical calculation using (54) and (32) produces the 
following typical results for the Dirichlet effective actions. For [3,3,3], W ~ 44.4 
and for [3,3,4], W ^ -427.25. 

10. Conclusion. 

The results of this paper arc strictly technical. We have achieved our aim of pre- 
senting calculable formulae for the functional determinants of minimally coupled 
scalar fields on the fundamental domains of finite refiection groups. The problem 
has devolved upon an evaluation of the derivative of the Barnes ^-function. 

We could also extend our previous results on the vacuum energies [1] to minimal 
coupling using the expressions for ({—n), (12), and ('{—n). This straightforward 
exercise will not be done here. 

The Cesaro-Fedorov formula for the number of symmetry planes of a regular 
solid proved in section 6 is one of a number of similar relations in higher dimensions 
derivable from expressions for the coefficients in the short-time expansion of the 
heat-kernel. The details will be presented elsewhere. 

The conformal transformations taking a fundamental domain into the upper 
half-plane are known and so the results here described should also be obtainable 
using standard conformal techniques. This will be recounted at another time. 
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Appendix. 

In this appendix wc first work out an expression for the derivative of the ^-function 
(8) at negative integers, ('{—n). For brevity we do not display the dependence of 
the Barnes (^-function on the d. 

Differentiation of (10) first of all leads to 



n r n ^ 

':'i-n)^J2(-^'yCr) 2Cd(2r - 2n, a) - Cd(2r - 2n, a) Yl k 

r—O L k=n—r+l 



k=n—r+l 
r—n—1 



+ 



(-1)" "f U-,.+iN.r-.:Ei)+ 



r=n+l 



fc=n+l 



(-1)'^ 2^ a^^^ ^Cd(2r-2n,a) 



(55) 



We substitute the integral form of the Barnes (^-function into the last term and 
find it as the s —2n limit of 2"(— I)"'??,! times 



2 / exp(— ar) cosh ctr — — — -— I ——p ^ (56) 

Jo \ ^0 (2r)! ; a (1 - exp(-d,r)) ^ ^ 

which equals 

u+n 2r 

T{s) (Cd(s, a-a) + Cd{s, a + a)) - 2 J] Tyry r(s + 2r)Cd(s + 2r, a). (57) 
The pole cancellation gives the condition 

Cd(-2n, a - a) + Cd(-2n, a + a) - 2Cd(-2n, a) = 



2./^2n\,, , ^ 2r2r-2n-l)! 

r=l ^ ^ r=n+l ^ ''' 

Extracting the finite remainder yields, after using (58), 

^ (^Cd(-2n, a-a) + C^(-2n, a + a) - 2C^(-2n, a)^ + 



2r-2n- 



(2n)! ^ 



2r 



(2n) 
2n 



2r 



(58) 



(^'(l + 2n- 2r) - -0(1 + 2n))Cd(2r - 2n, a) - Cd(2r - 2n, a) 
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(2r)! 

r=n+l ^ ' 



{lP{2r - 2n) + ^^(1 + 2n))N2r-2n + i?2r-2n 



(59) 



Multiplied by 2"(— l)"^n!, (59) must be substituted into (55) to yield a calcu- 
lable formula for C,'{—n). Doing so reveals that the remainder terms i?2r-2n cancel 
but, apart from this, there are no other simplifications apparent and we leave the 
analysis at this point. 

We next derive the result (44) starting from (38) and (39) whence 

e-a(0)+Ci,(0) _ r2(c?o)r2(c?o + l) . . 

r2(i) • ^ ^ 

It is necessary to use some properties of the multiple F-function. 
From (26) and (37) it is obvious that 

limF^(a) = -. (61) 

a^O a 

The other properties we need follow from the important recursion formula satisfied 
by the Barnes (^-function, 

C,r{s,a-\- di I d) - Cr(s,a | d) = -Cr-i(s,a | d'), (62) 

where d' stands for the set of degrees d with the di element omitted. 
If this equation is differentiated, it quickly results that, [2], 



Tr{a) Tr-i{a) 



(63) 



Tr{a + di) Pr-i(d') 
Setting a equal to zero in (63) and using (61) it follows that 

Vr{di)^Pr-i{d'). (64) 

For the group [g], we recall that the degrees are d = (g, 1). Then, choosing 
di = di = q and setting a = 1, we have from (63) 

r2(i) ri(i) 



F2(l + g) pi(l) 

which is clearly independent of q since the quantities on the right-hand side are 
calculated for the degree d' = (1). Further, from (64), it is likewise clear that F2(g) 
is independent of q. Therefore the quantity in (60), 

r2(g)r2(9 + i) _ pjii) 
m) ri(i)' 
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is independent of q. The actual value is 27r, agreeing with the particular cases (42) 
and (43). 

Incidentally, from the general formulae (2), (3) and (62) it easily follows that 
the [q] conformal ^-functions are related by 

Cf H«) - C'^Hs) = Ci(2., 1/2 I 1) = CRi2s, 1/2) (65) 

so that, in particular, 

Cf^'(0)-Cf^'(0) = -ln(2) 

for all [q\. 
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